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Abstract. This paper communicates recent results in theory of complex 
symmetric operators and shows, through two non-trivial examples, their potential 
usefulness in the study of Schrodinger operators. In particular, we propose a 
formula for computing the norm of a compact complex symmetric operator. This 
observation is applied to two concrete problems related to quantum mechanical 
systems. First, we give sharp estimates on the exponential decay of the resolvent 
and the single-particle density matrix for Schrodinger operators with spectral 
gaps. Second, we provide new ways of evaluating the resolvent norm for 
Schrodinger operators appearing in the complex scaling theory of resonances. 
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1. Introduction 

About half a century ago, Glazman laid the foundations for the theory of unbounded 
complex symmetric operators Since then, his fundamental ideas have been 

successfully tested on several classes of differential operators |7| 1201 Recently, 
two of the authors discovered an additional structure in the polar decomposition of 
a complex symmetric operator jl4j. For certain unbounded operators with compact 
resolvent, the refined polar decomposition leads to a new method for estimating the 
norm of their resolvent. In the present note, we exploit this idea in conjunction with 
the complex scaling method for Schrodinger operators. 

Although Quantum Mechanics is built on the theory of selfadjoint operators, 
one must often deal with non-selfadjoint operators. For instance, this is the case 
when appealing to the complex scaling technique. This method became a standard 
tool in the theory of Schrodinger operators and turned out to be the key to several 
problems such as: the absence of singular continuous spectrum the calculus 

of resonances and the convergence of time-dependent perturbation theory [23 , and 
asymptotic behavior of the eigenvectors jJUj. As the examples in this paper show, the 
complex scaling technique naturally leads to complex symmetric operators. 

Our first application deals with Schrodinger operators with a spectral gap. We 
provide sharp exponential decay estimates on the resolvent and the single-particle 
density matrix. Such estimates became increasingly important since it was realized 
that the localization of the single-particle density matrix provides the key to efficient 
numerical electronic structure algorithms for systems with large number of particles 
PSj- For ID periodic insulators, exact exponential decays can be derived from Kohn's 
analytic results [2]. In an attempt to generalize these results to higher dimensions, 
des Cloizeaux [2110] developed a method which can be regarded as the first application 
of the complex scaling idea. He proved the exponential decay of the single-particle 
density matrix for a class of 3D insulators. Relatively recently, we have seen a renewed 
interest in the subject and remarkable new exact results in dimensions higher than 
one ^ El El ESI- These results, however, are limited to periodic systems and some of 
them to the extreme tight-binding limit. In the present note, we treat the general case 
of gapped Schrodinger operators, which find applications, in addition to the periodic 
insulators, to amorphous insulators, molecular liquids, or large molecules. 

For periodic systems, it has long been known that the exponential decay of the 
resolvent, when the energy is in the gap and close to the gap edges, is proportional to 
square root of the distance from E to the gap edge. The theory of the effective mass 
\n\ provides a simple way of estimating this exponential decay constant. Relatively 
recently, it was proven that this qualitative behavior is present in any gapped systems 
PlEj- The present paper sharpens these previous estimates. The goal is to find 
the best quantitative estimate of the exponential decay constant, with the energy 
spectrum as the only input. 

In our second application, we show that the technique of estimating norms 
of complex symmetric operators extends to certain operators with non-compact 
resolvent, such as the complex scaled Hamiltonians from the problem of resonances. 

2. Complex symmetric operators 

This section is a brief account of relevant results (both old and new) about complex 
symmetric operators. For full details and examples the reader can consult |13[I14| . 
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We first consider bounded operators. Let Ti denote a separable complex Hilbert 
space which carries a conjugation C : H — > Ji, an antilinear operator satisfying the 
conditions = I and {Cf, Cg) ~ {g, /) for all f,g£H. A bounded operator T on H 
is called C -symmetric if T — CT*C. More generally, T is called complex symmetric 
if there exists a C such that T is C-symmetric. The terminology arises from the fact 
that T is complex symmetric if and only if it has a symmetric matrix representation 
with respect to some orthonormal basis 13' . 

Examples of bounded complex symmetric operators include normal operators, 
Hankel operators, finite Toeplitz matrices, Jordan model operators (the infinite 
dimensional analogs of Jordan blocks), the Volterra integration operator, and several 
other classes as well (see ^l^j). Some of these classes of operators appeared in direct 
applications like atomic collisions P and bound states problem in periodic chains [2^1 ■ 

The following simple factorization theorem (from |14j ) is the main ingredient in 
the proofs contained throughout this note. 

Theorem 2.1 \i4Jl IJ T : TL — > TL is a hounded C-symmetric operator, then there 
exists a conjugation J : Ti. — > Ti., which commutes with the spectral measure of 
\T\ = Vt*T, such that T = CJ\T\. 

The theorem above asserts the equivalence of the antilinear eigenvalue problems 

Tf = \Cf ^ \T\f = \Jf. (1) 

We may assume that A is real and positive, since we may multiply either of these 
equations by a suitable unimodular constant. This equivalence is important for the 
following reason. The norm of a compact operator T is equal to the largest eigenvalue 
of |T|. However, working with |T| instead of the original operator T introduces 
unwanted and potentially serious complications. Eq. provides the following 
convenient alternative: the norm of a compact C-symmetric operator T is equal to 
the largest positive solution A of the antilinear eigenvalue problem, T f — XCf. This 
is because the auxiliary conjugation J fixes an orthonormal basis of each spectral 
subspace of \T\. 

We now turn to the case of unbounded operators. A densely defined, closed 
operator T is C-symmetric if T C CT*C and C-selfadjoint if T = CT*C , i.e. their 
domains coincide and, on this common domain, they are equal. 

The study of unbounded complex symmetric operators was pioneered by Glazman 
|15l who established a complex symmetric parallel to von Neumann's theory of 
selfadjoint extensions of symmetric operators, although certain classes of C-symmetric 
operators had appeared earlier in von Neumann's work [34| . A renewed interest in 
Glazman's theory was sparked by its application to certain Dirac-type operators [7] 
and the realization that the closely related class of C -unitary operators is relevant to 
the study of complex scaling transformations in quantum mechanics [221 ■ Moreover, 
certain Sturm-Liouville operators with complex potentials can also be treated similarly 
| 20II27| . Further examples are furnished by Schrodinger operators —l^-\-v with complex 
potentials v (where C is simply complex conjugation) subject to appropriate boundary 
conditions One can also consider Schrodinger operators — A -I- w with real 

potentials v, but complex (non-selfadjoint) two point boundary conditions, in which 
case the conjugation C is slightly more involved. 

From the classical theory of selfadjoint operators, one knows that a symmetric 
operator has selfadjoint extensions if and only if its deficiency indices are the same. In 
contrast, every C-symmetric unbounded operator admits a C-selfadjoint extension T. 
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Unfortunately, not all unbounded C-selfadjoint operators possess a spectral resolution 
and a corresponding fine functional calculus. Nevertheless, if an unbounded C- 
selfadjoint operator has a compact resolvent, then a canonically associated antilinear 
eigenvalue problem always has a complete set of mutually orthogonal eigenfunctions: 

Theorem 2.2 \14\ If T : I?(r) — > H is an unbounded C-selfadjoint operator 
with compact resolvent [T — z)^^ for some complex number z, then there exists an 
orthonormal basis {un)^^i of Ti. consisting of solutions of the antilinear eigenvalue 
problem: 

(T ~ z)Un = XnCUn 

where {Xn)^=i is an increasing sequence of positive numbers tending to oo. 

This result is a direct consequence of the refined polar decomposition T = CJ\T\ for 
bounded C-symmetric operators described in Theorem 2.1. Our main technical tool 
in estimating the norms of resolvents of certain unbounded operators is contained in 
the following corollary: 

Corollary 2.3 \i4Jl If T is a densely- defined C-selfadjoint operator with compact 
resolvent [T — z)~^ for some complex number z, then 

\\iT-z)-^\\ = -^ (2) 
mf „ A„ 

where the A„ are the positive solutions to the antilinear eigenvalue problem: 

{T - z)Un = XnCUn. (3) 

Here \\ \\ denotes the operator norm: \\A\\ = supn^n^]^ {A(j), Acj)) . 

We also remark that the refined polar decomposition T = CJ\T\ applies, under 
certain circumstances, to unbounded C-selfadjoint operators. The following extension 
will play a central role in our second application: 

Theorem 2.4 \14Jl IfT : D{T) > Ti. is a densely defined C-selfadjoint operator with 

zero in its resolvent, then T — CJ\T\ where \T\ is a positive selfadjoint operator (in 
the von Neumann sense) satisfying T>{\T\) = T>(T) and J is a conjugation operator 
on Ti. which commutes with the spectral measure of \T\. Conversely, any operator of 
the form described above is C-selfadjoint. 

We close this section with a few remarks about the computation of the eigenvalues 
An = A„(|T|) and associated eigenfunctions via the antilinear problem or ((^J. If 
T is a complex C-symmetric operator, then the bilinear (as opposed to sesquilinear) 
form 

[Tx,y] = {Tx,Cy) = {\T\x,Jy) 

can detect these values via a min-max principle analogous to the corresponding 
procedure for selfadjoint operators (see |2Si XIII. 1]). More specifically, if Aq > Ai > 
X2 > ■ ■ ■ > are the eigenvalues of the compact positive operator \T\, arranged 
in decreasing order and repeated according to multiplicity, then the evenly indexed 
eigenvalues A2n are given by the following variational principle: 

X2n — min max ^[Tu,u]. (4) 

codimV— n -(iGViHiiH — 1 

The same principle applied to a rank one perturbation T -I- ^ (E) ^ (where ^ is any vector 
satisfying ||^|| > ||T||) yields the odd eigenvalues. The proof of Q), along with several 
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applications, can be found in We also remark that generalizes the recently 
discovered min-max principle for symmetric matrices with complex entries 

If one is only concerned with || T ||, then simply take n = in Q and note that 

||r||= max mTu,u\. 

new, ||u||=l 

The antilinear eigenproblem T f = AC/ and min-max principle 10} can be broken into 
real linear problems and also have the advantage that they do not involve the direct 
calculation of T*T . Moreover, one might potentially find numerical approximations 
to the A„ via procedures similar to those in i28| XIII. 1]. 

3. Exponential Decay of the Resolvent for Gapped Systems 

In this section, we consider the problem of finding sharp estimates on the exponential 
decay of the resolvent for Schrodinger operators with a gap in the spectrum. A short 
account on the subject has been already given in the Introduction. 

We now formulate the problem and the main result. Let — denote the Laplace 
operator with zero (Dirichlet) boundary conditions over a finite domain (with smooth 
boundary) Vl C R''; let w(x) be a scalar potential, which is V|,-relatively bounded, 
with relative bound less than one. Throughout this section, all potentials v are 
presumed to be bounded from below. By measuring the energy from the bottom 
of the potential, we can assume without loss of generality that w(x) > 0. We also 
include a magnetic field, described by a smooth vector potential A(x). The total 
Hamiltonian is 

Ha ■■ 'DiVD — > L^{n), Ha = -(V + iA)^ + w(x). 

The assumption on Ha is that its energy spectrum a consists of two parts, a C 
[0, E^] U [£'+, oo), which are separated by a gap G = E+ — i?_ > 0. We refer to the 
spectrum a± above/below the gap as as the upper/lower band. The corresponding 
spectral projectors are denoted by P±. 

Let E e {E-,E+) and Ge = [Ha - E)'^ be the resolvent. We are interested in 
the behavior of the kernel G'£;(x, y) for large separations |x — y|. Instead of looking 
directly at the pointwise behavior, we take the average 

Gi5(xi,X2) = ^ J dx J dyG'£;(x,y), 

x-xi|<e |y-X2|<e 

where is the volume of a sphere of radius e in R''. The main result of this section 
is stated below. 

Theorem 3.1 For q smaller than a critical value qc{E), there exists a constant Cq^E, 
independent of Q, such that: 

|Gb(xi,X2)| <C,,£;e"«l^i-^^l. (5) 

Cq^E is given by: 

■^'■^ mm\E±-E-q^\' l-q/F{q,E) ^' 



Norm estimates of complex symmetric operators applied to quantum systems 
with 



The critical value qc{E) is the positive solution of the equation q = F{q, E). 

Proof. If Xx denotes the characteristic function of the e ball centered at x (i.e. 
Xx(x') = 1 for |x' — x| < e and otherwise), then one can equivalently write 

(5e(xi,X2) = (j~^{x^^,{Ha - -B)"^Xx2)- 

Given a vector q e R'' = |q|) of arbitrary orientation and magnitude, let Uq denote 
the following bounded and invertible map 

[/q : ^ L-'in), [C/q/](x) = e'^-fi^), 

which leaves the domain of .^a unchanged. Let -ffq,A = Uc^Ha.U^^ be the family of 
scaled Hamiltonians. Explicitly, they are given by 

ifq,A : 2?(V|,) ^ ifq,A = ^?A + 2q(V + iA)-g2. (8) 

We notice that for q 7^ 0, these are non-selfadjoint and are not even complex symmetric 
operators (with respect to any natural conjugation). We also note that the identity 

C/q(HA - E)-^U-^ = (Hq,A - E)-^ 

holds for all q € R"^ (this happens only for finite O). If = R'', the identity holds as 
long as the continuum spectrum (which moves as q is increased) does not step over 
E. We denote 

7(g,i?)= sup ||(7Iq,A-i?)-i. (9) 

q|=9 

As the following lines show, the entire problem can be reduced to estimating j{q,E). 
Indeed, if (pi(x) = e~'^^^~^^^Xxii^) and (/?2(x) = e'^^^~^^^Xx2i^)j then 

|Gi,(xi,X2)| = coT^liifi, iH^,A - i?)-V2)|e-'»("^-"=) 

< a;7ie2«^7(g,£)e-'»(^^-^=). 

Choosing q parallel to xi — X2, we infer at this step that 

|G£(xi,X2)| < iV-^e^'''j{q,E)e-''\^'-^^\. 

Next we estimate j{q, E) using the results on C-symmetric operators presented 
in the previous section. As we already mentioned, H^ x is not complex symmetric. 
However, note that the operators H^^j^ and iJ-q.A are dual: * a ^ ^-q.A- 
Furthermore, if C denotes complex conjugation, then CH^ x = Hq^ _pjZ. We define 
the following block- matrix operator H and conjugation C on Z/^(Q) ® L'^{Q.): 

It is a simple task to check that H is C-selfadjoint: H* = CHC. Moreover, 

||(H - E)-^\\ = ||(/fq,A - E)-^ = ||(if-q,-A - ^)-^||. (10) 
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According to the previous section, the antihnear eigenvalue problem (with A„ > 0) 

(H - £;)0„ = A„C0„ (11) 
generates an orthonormal basis 0„ in L^(ri) © i^(ri) and 

ll(H-£;)-^ll-— (12) 

iiiin„ An 

If we write 4>n = fn ® 5n , the antilinear eigenvalue problem Eq. Hll|) , after we apply 
a complex conjugation on the second equation, is equivalent to 

{Ha, A — E)fn = Xudn ^-[^g-j 

{H-q.A — E)gn ~ A„/„. 
With q small, such that E + q^ lies in the spectral gap, the polar decomposition 

Hj^-E-q^ =S\HA-E-q\ 

holds, where S — P+ — P-. We take the scalar product of the first equation in Eq. (|13|l 
against the vector S'/„. Keeping only the real part of the result and solving for A„, 
we find 

, _ !(/„, \H^-E- q^\U) + 2jR(5/„,q(V + iA)U)\ 

mSfnr9n)\ ■ ^ ^ 

After elementary manipulations, the second term in the numerator can be rewritten 
as 

5R(5/„, q(V + iA)U) = 23fi(/„, P+[q(V + ^A)]P_/„). 
Moreover, denoting 

Bq = P+\Hjs^-E- qX^I'^HW + iA)]|FA -E- g2|-i/2p_, (15) 

one obtains 

|(/„,P+[q(V+*A)]P_/„)| < \\\Ba\\{fnAHA-E-q^\fn). 
Eq. (dH) implies 

\n>mm{\E±-E-q'\}{l-2\\Bq 



2| 1 /I oil D \\\ 



/n||||.9n|| 

Similarly, by taking the scalar product of the second equation of Eq. \V6\ against Sfjm 
one obtains: 

A„ > min{ \E^-E-q^\ }(1 - 2||i?q||) 

The sum of the last two equations yields 

A„ > min{ \E±-E-q^\ }(1 - 2||Bq||). (16) 

It remains to evaluate ||Bq||- Since the potential is positive and we work with 
zero boundary conditions, the following inequality between quadratic forms 

q2(-(V + iAf +v + a)> [q(V + iA)f , V a > 0, 

holds true. Consequently 

||[q(V + iA)]|iJA + ari/2|| <g, V a > 0. 
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We can then insert \Hp^ + a\ ^/^|i?A + a|^^^ after q(V + iA.) in Eq. |{T^ and use the 
above estimate and the spectral theorem to get: 



Finally one can pass to the limit a ^ 0. Eq. H12|l . together with Eqs. H16() and (|17|l 
provide the desired estimate of ^{q,E).U\ 

We have the following remarks. The reason we considered only finite domains is 
primarily because the abstract machinery presented in the previous section applies only 
to complex symmetric operators with compact resolvent. An extension to operators 
with non-compact resolvent is given in the next Section. However, since the constants 
in Theorem 3.1 are independent of the domain fi, the estimates remain valid when 
R'', whenever this limit is well defined. For E close to the gap edges, a 
brute perturbation theory on (^fq,A — E)~^ leads to an exponential decay constant 
proportional to \E± — E\. The correct behavior of qc{E) near the band edges is 
\E±-E\^/^, as it was first shown in Ref. 3_ and later reviewed in Ref. US|- Note that 
our estimate of qc{E) does have the correct behavior near the gap edges. 

Corollary 3.2 Consider the lower band completely filled. Then the single-particle 
density matrix (i.e. the projector onto the occupied states P^) decays exponentially, 
with a rate q satisfying 

q>-^- (18) 



Proof. Again, we look at the average 

P_(xi,X2) = CJ72(Xxi,-P-Xx2), 

which has the following representation: 

P„(xi,X2) = ^^GB(Xi,X2)d£;, 

where F is a contour in the complex energy plane, surrounding the lower band. The 
estimates given in the preceding Theorem trivially extend to the case of complex 
energies: 

|G£;+.c(xi,X2)| < C,,Be-*l"i-"^l, V q < qc{E). 

Given that F can be deformed so as to intersect the real axis at any point in (-E-, i?+), 
we need to find the energy where qc{E) is maximum. We have 



= (E.-E-qmE-E^+q^) < G 



with equality at energy 

E^^±±^-^. (20) 
2 16E^ ^ ' 

The lower bound of Eq. H18|) is valid as long as E is in the gap.D 

The lower bound on q given in Eq. H18|) can be calculated entirely from the energy 
spectrum. What is the best lower bound on the exponential decay constant that one 
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Figure 1. The exponential decay constant of the single-particle density matrix 
as a function of G/W for the Kronig- Penney insulator. The continuous line 
represents an exact calculation and the dashed line represents the estimate given 
in Eq. 1181 . The inset shows the relative difference between the two. 

can get by only using the information contained in the energy spectrum? We illustrate 
by an example that Eq. H18|l comes close to such an optimal estimate. In Fig. 1, we 
consider a comparison between our estimate Eq. (|18|l and the exact value of q, for a 
one dimensional insulator described by the Kronig-Penney model |23| : 

H = -dl + Wo ^ 5{x - n), Wo > 0, 

n 

with the first band completely filled. The Dirac-dclta potential is not relatively 
bounded so the above example is not formally covered by Theorem 3.1. However, 
this potential is a relatively bounded form perturbation and our arguments can be 
extended to this larger class of potentials. By varying the strength of the potential wq, 
we sweep from a weak to strong insulating regime, which we quantified by the ratio 
between the gap G and the width of the valence band W . One can see that even in the 
extreme insulating regime (typically G/W < 5), Eq. I|18(l estimates the exponential 
decay to within a 15% error. Notice that for G/W < 5, the error is less than 5%. 

4. Norm Estimates on Resolvents Near Resonances 

In this section, we extend the results concerning norm estimates to not necessarily 
compact resolvents of unbounded complex symmetric operators. At the same time, 
we apply this technique to the problem of locating the resonances of a specific class 
of Hamiltonians. 

Wc first formulate the problem in precise terms. Let 



H : V{V^) — > L^{R''}, H = -V^ + t;(x) 
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be a Hamiltonian with w(x) a dilation analytic potential in a finite strip < /g s-^d 
V^-relatively compact. We consider the usual dilation operation: 

[{/(e)V'](x) ^e'^^/Vle^x) 

and define the analytic family (of type A) of operators: 

He = U{0)HU{9)-^ = -e^^^y^ + t;(e^x), 

where 9 runs in the finite strip < la. As a function of 6, it is well known that 
001211: 

(a) the discrete spectrum ad remains invariant, 

(b) the essential spectrum aess rotates down by an angle —23^?, 

(c) as the continuum rotates, it uncovers additional discrete spectrum (the 
resonances) . 

In many practical situations, it is desired not only to locate the resonances but 
also to know how they move under different perturbations |31 1^ 123 ISH ■ Here we are 
concerned with the second problem, where norm estimates on the resolvent {z~ Hg)~^ 
for z near the resonances become especially important either for probing the stability 
of the spectrum or for building perturbation series. 

The Hamiltonians Hg are C-selfadjoint relative to the complex conjugation 
Cf = f. The question that we want to answer is if one can provide an exact norm 
estimate of {z — Hg)^^ for z near a resonance, using the theory of complex symmetric 
operators. The answer is contained in the following theorem: 

Theorem 4.1 Let jw{:x.) represent the change in u(x) and H{'-f) = H + 71/7 denote 
the perturbed Hamiltonian. We assume that both vg{x.) = w(e^x) and Wg{x.) = w{e^x) 
are 'S/'^ -relatively bounded for < Iq, with bound less than one. For z close to a 
resonance zq of H and 7 sufficiently small, the following are true: 

(l) CJess{\Hg{-l) ~ Z\) = [d(z,0),oo). 

(li) Ud{\Hg{-i)- z\)n[Q,d{z,9)]^%. 

(Hi) Xn G ad{\Hg{'j) — z\) if and only if there exists ijjn G Viy^) such that: 

{Hg{j) - z)%lJn = XnClpn- 

Moreover 

\\{Hg{^)^z)-^\\^-±—. 



Above, d{z,6) = |zsin(230 — a)\ denotes the distance from z to aess{He) (where 
z = |z|e-*";. 

We require the following lemma: 

Lemma 4.2 Let A and B be two closed operators such that'D{A) C 'D{B) and B\A\~^ 
is compact. Let A + B be the closed sum on ^{A). Then a^sslA + B\ — aess\A\. 

Proof. Let C ^ \A + B\^, defined on \A\-^n. We show that (C - C^)~^ is a 
meromorphic operator valued function on C G C\[<t(|A|) U(7(— |^|)]. This follows from 
the identity 

(c - er' = {\A\ + c)-Mi + A^(c)]-'(i^i - c)-' (21) 
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where 

Ar(C) = {\A\ - C,)-'^[A*B + B*A + B*B]{\A\ + C)"^ 

is an analytic family of compact operators on ^ G C\[cr(|A|) U cr(— |A|)].n 

Proof of Theorem 4.1. i) Taking A = — e~^^V^ — z and B = vg + -ywg, it follows 
that the essential spectrum of \Hg{'-f) — z\ is contained in cr(| — e~^^V^ — z\), which is 
[d(^,^),oo). 

ii) Wc need to show that |-ff6i(7) — z\ has spectrum below d{z, 0). Let f/'o be the 
eigenvector corresponding to the resonance, HgipQ = zqiJjq. Remark that 

(Heij) - ^)Vo = {zo - z)[l + jweiHe - z)-^]ipo 

and consequently 

\\\He{j) - zlM < \zo-z\{l + j\\we{He-z)-'\\). 

With our assumptions, there exists < a < 1 and 6 > such that IjweV'll ^ 
a||V^VII + Hi^W ^^'^ similarly for ve, for any ^ e ©(V^). A relatively elementary 
manipulation then yields: 

The conclusion is that we can make |||-ff0(7) — 2;|V'o|| arbitrarily small, in particular, 
smaller than d{z,9), by taking the limit z ^ zq and 7 — > 0. Consequently, 
ini a{\Hg{j) — z\) < d{z,9) for 7 small enough and z close enough to the resonance. 

Hi) Since the operator H${'-f) — 2; is C-selfadjoint, it admits the decomposition 
stated by Theorem 2.4: 

He{7)-z = CJ\He{7)-z\, 

where the second conjugation J commutes, in the strong sense, with the sclfadjoint 
operator \Hg{'y) — z\. In particular J leaves invariant the spectral subspaces of 
IHei'f) — z\. Thus if A„ belongs to the discrete spectrum of |i?0(7) — z\, then the 
vector space consisting of the eigenvectors e ©(V^): 

|-ffe(7) - Z\^n = K^n 

is left invariant by J. Thus, either ^„ = —J(pn or (ji^ = ^„ + J(^„ provide a new 
eigenvector tpn satisfying J'ipn = V'n- Therefore, 

C{He{j) - z)i>n = J\He{l) - z\i>n = \He{l) - ^|Vn = A„Vn- 
This proves the last assertion in the statement. □ 

5. Conclusions 

The goal of this paper was to introduce recent results in theory of complex symmetric 
operators and to show, through two non-trivial examples, their potential usefulness 
in the study of Schrodinger operators. It was shown that, for complex selfadjoint 
operators with compact resolvent, a certain antilincar eigenvalue problem can be 
used in a similar way the spectral decomposition for self-adjoint operators is used 
to estimate the norm of the resolvent. We combined this observation with the 
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complex scaling technique, to obtain sharp lower bounds on the exponential decay 
of the resolvent and density matrix for Schrodinger operators with a gap in the energy 
spectrum. We also shown that these techniques can be applied to complex selfadjoint 
operators with non-compact resolvent, in particular, to Schrodinger appearing in the 
complex scaling theory of resonances. 

In concrete applications, we believe that one can go beyond the approximations 
made in this paper and use either the newly developed min-max principle for complex- 
symmetric operators E]; or direct, numerical integration of the anti- linear 
eigenvalue problem to obtain sharper, problem specific estimates of the exponential 
decay constant or to evaluate the resolvent norm in our second application. 
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